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Abstract 

In this paper, a new filtering method is presented for simultaneous noise reduction and enhance- 
ment of signals using a fractal scalar conservation law which is simply the forward heat equation 
modified by a fractional anti-diffusive term of lower order This kind of equation has been first intro- 
duced by physicists to describe morphodynamics of sand dunes. To evaluate the performance of this 
new filter, we perform a number of numerical tests on various signals. Numerical simulations are 
based on finite difference schemes or Fast Fourier Transform. We used two well-known measuring 
metrics in signal processing for the comparison. The results indicate that the proposed method out- 
performs the well-known Savitzky-Golay filter in signal denoising. Interesting multi-scale properties 
w.r.t. signal frequencies are exhibited allowing to control both denoising and contrast enhancement. 
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1 Introduction 

Filtering is a process that removes some unwanted component from a signal. It is a very important task 
in signal processing, data analysis and communication systems. Many techniques have been proposed 
for this purpose. For instance, we can use simple averaging filters such as a Gaussian filter. It is well- 
known that this filter can be realized by solving the heat equation. Other more general partial differential 
equations (PDF) have been used, with non-linear anisotropic diffusion ||6l |22l, non-Unear fractional 
diffusion ||4l[15l or fractional time-derivative Q. It has been proved that PDFs are suitable in signal 
denoising EOll . The denoising applications have to take into account two points. First, we want to obtain 
a clean and readily observable signal (improving signal-to-noise ratio SNR) and secondly, preserve the 
original shape characteristics (maxima, minima...) of the signal. This task is complex because it is very 
important that the denoising has no blurring effect on the images and does not change the location of 
image edges, fSj. For some applications, it is interesting to amplify some features of the signal such 
as relative maxima or minima in order to enhance its contrast. Usually, these features are flattened by 
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the denoising methods based on averaging techniques. Among the denoising methods which preserve 
characteristics of the initial signal is the Savitzky-Golay filter, with which we will compare our new 
method. 

A basic and crude idea to enhance the contrast could be to solve the backward heat equation for a 
few time steps. Of course, this is an ill-posed PDE and we do not advocate this unsafe method but it 
illustrates the fact that enhancement and denoising are antagonistic operations. Our method is based on a 
linear PDE, with two antagonistic terms : a usual diffusion and a nonlocal fractional anti-diffusive term 
of lower order. The diffusion is used to reduce the noise whereas the nonlocal anti-diffusion is used to 
enhance the contrast. Let us emphasize that we perform at the same time noise reduction and contrast 
enhancement. Our method is based on the Cauchy problem of the following PDE : 



{dtu{t,x) - adl^u{t,x) + hlx[u{t, .)]{x) = t G (0,r),x G 
1 n(0, x) = uq{x) X E M, 



(1) 

where T is any given positive time, -uo G L^(M),a, 6 are positive constants and Xa is a fractional operator 
defined as follows via the Fourier transform: for any Schwartz function tp G SiM) 

ZxM{x):=-F-\\.\^H^)){x) (2) 

where < A < 2 and F denote the Fourier transform defined by: for all ^ G M 

Let us note that for A G]1, 2[, we have an explicit nonlocal formula (see proposition [5]) 

I,Mx)=a, [ ^^^de (3) 
Jr I?I 

where ax is a suitable constant. 

Alternatively, we can also give a slightly different definition, inspired by fractional calculus (see remark 
below) : 

^"{x-O 



ZxW\{x)=axj^ "|^|A„i (4) 







Remark 1. For causal signals (i.e ip{x) = 0/or x < 0) we have 



1 d^""^ //_ 



r(2 - A) 7o ^ d^"" ^ d^"" 

which is exactly the Riemann Liouville definition of the fractional derivative 

Remark 2. Our model is closely related to a nonlocal conservation law first introduced to describe 
the morphodynamics of sand dunes and ripples sheared by a fluid flow. Namely, Fowler (/ [72] 173] 1741/ ) 
introduced the following equation 

dtu{t, x) + {t, x) + I[u{t, .)]{x) - dl^u{t, x) = 0, (6) 
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where u = u{t, x) represents the dune height and X is a nonlocal operator defined as follows: for any 
Schwartz function ip G and any x € R, 



Equation ^ is valid for a river flow over a erodible bottom u(t, x) with slow variation. The nonlocal 
term appears after a subtle modeling of the basal shear stress. See /[7] ^ for theoretical results on this 
equation. 

This nonlocal term appears also in the work of Kouakou & Lagree / ITTl 17^ . The operator X[u] is a 
weighted mean of second derivatives ofu with the bad sign and has therefore an anti-diffusive effect and 
create instabilities which are controlled by the diffusive operator —d"^^. We find again this phenomenon 
for the equation ([T]l. 

The solution of the Unear PDE ^ is obtained by convolution with the kernel KofZx — d'^^. There- 
after, this kernel will be our filter for denoising and enhancement of signals. The analysis of this kernel 
shows that the low frequencies are preserved, the medium frequencies are amplified and the high fre- 
quencies are eliminated. It is cleai" that this kernel depends on the parameters a, b, X and that the choice 
of these coefficients will determine the quality of the noise reduction and of the enhancement. In this 
paper, we discuss the choice of these parameters. 

To evaluate the performance of our filter, we discretize the PDE by two methods : the fast Fourier trans- 
form and the finite difference method. Numerical studies of nonlocal equations ai^e rather scaixe: among 
them, we mention lITTI which proves the convergence of a finite volume method to approximate the so- 
lutions of a fractal scalar conservation law, that is to say a conservation law regularized by a diffusive 
fractional power of the Laplace operator and Q which analyzes the stability of finite difference schemes 
for the solution of In these works, the discretization of the fractal operator is performed via an 
integral formula for Ix similar to the ones appearing in |[T]|9l- Our method is then compared with the 
Savitzky-Golay filter. Results show that the PDE ([T]) is relevant and effective for denoising with preser- 
vation or enhancement of features of the signal. 

The remaining of this paper is organized as follows: in the next section, we explicit the solution 
of problem ([T|l and we give some properties of the kernel. In section [3l we present explicit numerical 
schemes which approximate the fractal conservation law ([T|) and we give some numerical simulations. 
In this section, we also briefly present the Savitzky-Golay filter and compare it with our filter. Other 
numerical simulations based on the EFT are given in section HI We also discuss the choice of parameters 
a, b, X and we highlight both the ability of denoising and contrast enhancement of our model. Section |5] 
is devoted to the performance evaluation and metrics. Section |6] gives conclusions about this study. 

2 Theoretical study of the PDEs 

In this section, we verify that ([T]) is well-posed and in subsection 12.21 we analyze the properties of the 
kernel KofXx- dl^, for < A < 2. 

2.1 Well-posedness of the problem 

Using the Eourier transform, we see that any solution to ([T]) satisfies the formula dSjl. 




(V) 
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Proposition 1. LetT > and uq G L'^{R). The function u G L°°{{0,T); L^{R)) is a solution of ^ if 
for any t G (0, T): 

u{t,x) = K{t,-)*uo{x) (8) 
where K{t, x) = T^^ (e"*'^^'^) (x) with ^{5^ := A^'K^a^ — b \^\^ is the kernel of the operator Xx — d^^. 

Proposition 2. LetT > 0, uq G L^(IR). Then, the function 

u : t G (0,r] ^ K{t,-)*ua 

is well-defined and belongs to C([0, T]; Lp'iM)): u is extended at t = by the value u(0, .) = uq. 

Proof. We have 

K{t,x)= [ e^'^<h{i)di, 

where h{^) = e'^'^^?). For < A < 2, it is easy to verify that h belongs to ^j^ere W^'^{R) := 

G 1^, G L^I^)} therefore we have that K(t, .) G L^(M). Hence, Vt > 0, K{t, •) * uq is 

in L2(M). 

Let us prove the strong continuity ie 

lim K (t, .) * uo = Uq in (M) . 

t-s-o 

By Plancherel's formula, 

I •) *'Uo - nol 1^2(18) = \ \T{K{t, ■) * Uo) - Tuo\\l2^f.) 

= \\e-'^'J^uo - -^^xo||i2(K) = / le-*'^ - Ip I-FtxoI'. (9) 

Since the function \e~''^^ ~ |-^^oP converges pointwise to on M, as t — )• and as min^ is finite 
then, by the dominated convergence theorem, the last term of Q tends to as t 0. 

■ 

Remark 3 (Regularity of the solution). It is easy to see that u G C°°((0, T] x M) because K is smooth. 
The smoothness of K is an immediate consequence of the theorem of derivation under the integral sign 
applied to the definition of K by Fourier transform. To obtain the regularity att = 0, we have to suppose 
that the initial condition uq is in C°°(M) and satisfies for all /c G N, G L^(M). The proof is similar 

(k) 

with replacing uq by Uq . 

Remark 4. From formula (H), it is straightforward that 

/ = J"(A'(i,-)) (0) • / uo{x)dx= / uo{x)dx. (10) 

This is a "mass conservation " property. 



4 



2.2 Study of the kernel 

In this subsection, we give some properties on the kernel Kofix — d"^^. 
Proposition 3. The kernel K has a non-zero negative part. 
Proof. Let us assume that K is nonnegative, then 

|e-#(0| < ||^"i(e-*'^)||^.(^)= / \K{t,.)\ 

Jr 

= [ Kit,.)=TiK{t,.))iO) = e-'^^'^ = 1 
Jr 

for all C G M. But we have for < |C| < {j^)^ , le-^'^^?)] = e-*(^'^'<'-^l«l^) > 1 , this yields a 
contradiction. ■ 

The main consequence of the non-positivity of K is the failure of the maximum principle for the 
equation ([1]) IH. Thereby, u{x, t) is not forced to remain in the interval [min(iio), max(iio)]. The signal 
can then be amplified when it is necessary. 

Using the proof of the previous proposition, we can deduce that the enhancement of the frequencies will 
be feasible only for the low/medium frequencies. Of course, this amplification will depend on the choice 
of the parameters a, b, X and of the time t. 




■-8 -6 -4 -2 2 4 6 8 5 10 15 20 25 



(a) The kernel of Xx for t = 0.1 (red) and t = 0.5 s (blue) (b) Behaviour of tp 

Figure 1: a = 0.5, 5 = 2 and A = 1.5 

We expose in figure [T] the evolution of K{t, .) for different times and the behaviour of if; for a = 
0.5,5 = 2, A = 1.5. 

3 Finite difference schemes 

In this section, we present a finite difference numerical scheme to directly approximate the solution to 
([B for any A G]0,2[. 
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3.1 Integral representation of Xa 

To approximate the nonlocal term, it is useful to give an integral formula for X\: 
Proposition 4 (||T0l[I3). //A g]0, 2[, we have for all tp G cS(M), a// x G M and all r > 0, 

/ r 99(x + 2:)-v?(a;)-99'(x)2: f (^(x + z) - v?(x) ^ 

x.bi(x) = (^1^ + — — 

l+A \ 

where C\ = — Ta"^ — '^'^'^ ''^ denote the Euler function. 
1) IfX £]0,l[then 

XAb](a;)=CA / ^-^^ dz. 

2j //A G]l,2[?/te?i 

r 1/ N ^ /■ 99(X + Z) - (/?(x) -(^'(x)2: 

2^AM(a;) = Cx / ^ta dz. 

From this proposition, we deduce the following useful result : 
Proposition 5. For all tp G 5(M) and all x we have for A G] 1, 2[ 

f p (x ~\~ z\ f ip ix — 

X\[(p\{x) = ax I |^_. dz = ax / , dz 

JR 1^1 JR 1^1 

where ax~ ^ 



A(A-l)- 

Proof. The regularity of p ensures the validity of the following computations . Since 

f 2 



(p{x + z) - p{x) - LP {x)z = I (^p (x + y) - ip (x)^ dy 

p {x + tz) — p (x)^ z dt, 



the last equality arises from the change of variable y = tz. Then, using Fubini's Theorem, we have 

Xx{x) = Cx {p {x + tz) - p (x)) ^ dtdz 

JrJo 1^1 

( / + tz) - p'{x)) dz ) dt 

) \Jr / 



y=tz 



Cx I { I {p'{x + y)- 9^' (2;))^ dy ) t^-^ dt 



Cx I t^-Ut I ipix + y)-p{x))j-^dy 



\y\' 

\y\' 



z=sy 



\ I lr|A-l 
JR \Z\ 

Cx f ip"{x + z) 



3.2 The numerical scheme 

The spatial discretization is given by a set of points Xi;i = 1, ...,N and the discretization in time is 
represented by a sequence of times = < ... < < ... < T. For the sake of simplicity we will 
assume constant step size 6x and 6t in space and time, respectively. The discrete solution at a point will 
be represented by w u{t'^,Xi). In this section, we will present the behaviour of explicit numerical 
scheme which directly approximate the solution to ([T|). 

We discretize all terms of the equation using an explicit method. We consider for any n,j G N 

' St ' - « J — + bTs.K]j = (14) 

where Isx is a discretization of the nonlocal term Ix- This scheme is explicit because the values of the 
solution at time are obtained directly from the (known) values at time t". For the Laplacian term, 
we use a standard centered finite difference approximation of second order. To discretize the fractal 
operator Ix, we consider the formulation (|4]), which for A 2[ is a causal variant of ([T3] ). Next, we 
use a basic quadrature rule to approximate the integral and we use a finite difference approximation of 
the derivative: 

+ 00 

1=1 

Note that we have absorbed the constant ax in b. 

Remark 5. The practical implementation of the schemes requires to make truncations. The main trun- 
cation concerns the integral operator for the nonlocal term Xx- We replace J^°° with and the finite 
difference approximation becomes: 

A 

XsxWlj = 6x~^ ^/^~'^(-yj„i+i - 2vj^i (16) 
1=1 

with L = A5x. Usually, we consider that functions are either compactly supported or constant near 
— oo and +oo, so it is legitimate to consider a finite sum for the discretization of Xx- However the 
truncation parameter A has to be chosen judiciously. Indeed, when < A < 1 , the term in the 
discretization (1161 ) is bigger and bigger for increasing I, hence A ought to be big enough. In contrast, 
whenever 1 < A < 2, /^"^ is negligible for large values and is important only for small values of I. Thus 
in this case, it is judicious to take A small enough. We see again here that the non local effect is all the 
more important than A is small. These memory effects strongly depending on A have been reported in 
/[S]/. To take this behaviour into account, we set the truncation parameter A := max(100, Of course 
100 is an arbitrary limit to avoid rounding effects for small A. 

Since the equation is implemented using an explicit method, this imposes a CFL condition on the 
time and space steps which ensures the numerical stability that is to say that the difference between the 
approximate solution and the exact solution remains bounded when n — +oo for 6x, 6t given. The 
stability analysis of the nonlocal scheme (fT4l ) is done in [3J. This requires a careful study, owing to 
the anti-diffusive behaviour of the nonlocal term. Indeed, the analysis of the PDF ([T]) shows that for 
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the continuous problem, low to medium frequencies are amplified by the nonlocal term, IS. We come 
back to this in section \4\ Therefore, the standard Von Neumann definition of stability must be slightly 
modified. Denoting 9 = k 6x ^ {0, 2tt) and = exp {ij9), instead of imposing that the amplification 

factor G{9) = must fulfill |G(0)| < 1 for any frequency 9, we only impose < 1 for ^ above 

a given threshold. We obtain the two following conditions on the time and space steps: 

2a6t , x , b5t 

The first stability condition ([TTl) forces the mesh-size 6x to be small enough in order that the diffusion 
term should dominate the nonlocal anti-diffusive term for high frequencies, whereas the second stability 
condition (ITSl) looks like an usual CFL condition for explicit schemes and forces the time-step to be small 
enough. For more details, we refer the reader to ISj. 



3.3 Numerical results 

In the following numerical tests, we have to take care to choice 6t and 6x accuratly following the con- 
ditions (fTTl) and (fTSl ). Thus, the time and space steps depend on the choice of a, b and A. We begin 
by consider an electrocardiogram (ECG) signal. Figure |2] illustrates a comparison between the PDF ([D 
implemented using a finite difference method (FDM) and the Savitzky-Golay smoothing filter. 
The Savitzky-Golay smoothing filter also called digital smoothing polynomial filter or least-squares 
smoothing filter was first described in 1964 by Abraham Savitzky and Marcel J. F. Golay. |[24l . Fs- 
sentially, the Savitzky-Golay method performs a local polynomial regression on a series of values that 
is to say it replaces each value of the series with a new value which is obtained from a polynomial fit to 
neighboring points. The algorithm is based on the following equation 



n+l 
2 



1=1 



where n is the number of data points, are (positive) constants, v is the noisy signal and u defines the 
filtered signal. The main advantage of this approach is the preservation of features of the signal such 
as relative maxima, minima and width. Indeed, usually these characteristics are 'flattened' by classical 
averaging techniques such as moving averages. Figure |2] shows three plots: The noisy FCG signal, 
the smoothed signal (red) using the PDF filter ([T]) superimposed with the noiseless signal (blue) and 
the smoothed signal (red) using the Savitzky-Golay filter superimposed with the noiseless signal (blue). 
As we can remark, for these parameters, the low/medium frequencies are not amplified but the relative 
maxima and minima are preserved. The denoising seems con^ect and the filtering with our numerical 
scheme and with Savitzky-Golay seems similar, see better. We evaluate our approach and we compare it 
with Savitzky-Golay in section [5j 

Figure[3]illustrates two filtered signals of different types. We start from the simplest possible example 
of nonlocal filter denoising applied to a signal u{x) consisting of a piecewise constant step function 
v{x) = — 1 for < X < 1 and v{x) = 1 for 1 < a; < 2 corrupted by additive Gaussian white noise n{x) 
with standard deviation u = 0.4: 

u{x) = v{x) + n{x). (19) 
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Noisy ECG 
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Noiseless signal vs. filtered signal using FDM 
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Noiseless signal vs. Savitzky-Golay filter 
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Figure 2: Top: a noisy electrocardiogram u; Middle: in red the noiseless signal, in blue the filtered 
signal with finite difference scheme (fT4l) with a = 4, 6 = 0.5, A = 1.7, 6x = 1 and 6t = 0.1; Bottom : 
Noiseless signal (blue) vs. filtered signal using Savitzky-Golay filter. 



The result of the denoising is illustrated in figure [3(a)| As we can notice, the noise is very well eliminated 
and we find again our original signal, the shape of the signal has been preserved. The result is better than 
filtering Savitzky-Golay. Moreover, unlike Discrete Fourier Transform where the contrast is low in the 



neighbourhood of the discontinuity (see figure [5(a)| ), the dicontinuity/shock is conserved. Therefore, the 
finite difference method is more suitable for this kind of signal. Figure [3(b)] shows another example of 
good denoising. 

4 Numerical results using Discrete Fourier Transform 

In this section, we first discuss the choice of parameters a,b, X and next, we give some numerical results 
based on fast Fourier transform (FFT) which is an efficient algorithm to compute the discrete Fourier 
transform. 

4.1 Choice of parameters 

We fix T = 1. So, we can rewrite the kernel in Fourier variable as follows: 



where -0(0 = 47r^a^^ — &|^|'^. We draw in figured the behaviour of the kernel K^f^ for a, A fixed 
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Noiseless signal vs. Noisy signal 



Noiseless signal vs. filtered signal using FDM 



Noiseless signal vs. Savitzky-Golay filter 



iijju iukiiiiii 



Noiseless signal 



Filtered signal using FDM 



(a) a = 3.5, b = 0.2, A = 0.1, St = lO"*^ and (b) a = 3.5, b = 0.3, A = 1.1, St = 0.1 and Sx = 1. 
5x = 0.001. 

Figure 3: (a) Top: Noisy signal (blue) vs. noiseless signal (red); Middle: Noiseless signal (blue) vs. 
filtered signal using FDM (red); Bottom: Noiseless signal (blue) vs. filtered signal using Savitzky-Golay 
filter (red), (b) Top: Noisy signal; Middle: Noiseless signal ; Bottom: filtered signal using FDM . 
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(a) 47r^a = 0.01 and A = 1.5 (b) 47r^a = 0.01, 6 = 0.05 and A = 1.5 

Figure 4: Evolution of the spatial Fourier transform of the kernel K for different parameters. 



and K^JM 



and for different values of b. As K^^^ reaches it maximum at '■= (g^)^ ^ 

e'^^ ""^stv'^J ~ ^\~^\ we deduce that K^^{£^m) > 1 if and only if < A < 2. Therefore, whatever the 
choice of A g]0, 2[ , we always have an amplification of medium frequencies. Of course, the magnitude 



of amplification depends on a, 5, A. Another natural frequency is 



which represents the 



neutral frequency and satisfies K^^{^i) = 1 and K^^{^i) < 1 for ^ > ^i. Therefore is the threshold 
above which dampening will occur. 

Now, we wish to control both the denoising and contrast enhancement. For that, we follow a simple 
strategy: 



1. In a first step, we wish to control the frequency range to amplify. This one can be controlled by 
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the following ratio: ^ = (f)^^^ ^- Indeed, this ratio is decreasing w.r.t. parameter A and shrinks 
from +00 to when A varies from to 2. Hence, to choose a given amplified frequency range. 



1 

it is enough to fix parameter A and the ratio ^ (which in turn determines ^1 = (4;;^) ). 

2. In this step, we want to monitor the denoising, which will start for frequencies above the neutral 
one ^1- One can easily check the following equality: 



a 



where a is a any positive constant. Therefore, the greater a is, the smaller A'^j(a^i) will be, 
which means that the dampening rate will increase. Hence parameter a monitors the denoising 
intensity, this is coherent because it controls the Laplacian term. 

3. To finish with, we wish to monitor the value M := -f^afe(^Af) in order to control the contrast 
enhancement. Indeed, the higher M is, the better the constrast will be amplified. Parameters A 
and ^ being fixed by item 1 it is enough to adjust both parameters b and a while keeping the ratio 
^ within some bounds in order to preserve the amplified frequency range. The expression of M 
shows that to have a good amplification of medium frequencies, it is necessary to increase b as 
well as a. This is quite natural since the coefficient b controls the anti-diffusive term and has an 
opposite behaviour to the Laplacian term which tends to flatten the signal. In our method, contrast 
and denoising are no more antagonistic. We come back on the tuning of parameter b in section [S] 

4.2 Denoising 

In this subsection, we aie only going to highlight the ability of noise reduction of our nonlocal equation 

To begin with, we consider an ECG signal to illustrate the denoising using the FFT to solve the fractal 
equation ([T|l. The result is given in the figure |5(b)[ For this choice of parameters, we note that the 



Noisy signal vs. filtered signal 
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Noiseless signal vs. filtered signal 
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(a) 47r^a = 0.01, b = 0.03 and A = 0.5 



Noisy signal vs. filtered signal 
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Noiseless signal vs. filtered signal 
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(b) ATT^a = 0.0051, b = 0.015 and A = 1.7 



Figure 5: Top: Noisy signal (blue) vs. filtered signal (red). Bottom: noiseless signal (blue) vs. filtered 
signal (red). 

denoising is suitable and that the relative maxima are preserved. However, the relative minima are not 
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completely preserved, in fig. |5(b)| bottom. At last, Figure 10(a) [ describes the denoising of a sinusoidal 
signal. We can see the performance of the denoising of nonlocal approach discretized with FFT, which 
works very well for this type of periodic signal. 



4.3 Denoising & Enhancement of signals 

In this part, we are going to highlight both the ability of noise reduction and contrast enhancement. Of 
course, to emphasize the amplification, we take into account the remarks done in subsection 14.11 

We start from a simple example of simulatenous denoising and enhancement applied to a strongly 
attenuated sinusoidal signal highly corrupted by a random noise. In figure [6] middle, we plot the original 
noiseless signal, the same signal amplified x 50 and the filtered signal, performed with our non local FFT 
method. As we can notice in figure [6l the noise is eliminated and the contrast is well amplified. 



Noiseless signal vs. noisy signal 
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Noiseless signal (with or without amplification) vs. filtered signal using FFT 
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Figure 6: Top: Noisy signal (blue) vs. noiseless signal v{x) = 0.02cos(x) (red); Middle: Noiseless 
signal (blue) vs. filtered signal (red) using FFT with An'^a = 0.2, b = 1.25, A = 0.5 superimposed with 
amplified noiseless signal cos(x) (green); Bottom: Noiseless signal (blue) vs. filtered signal (red) using 
Savitzky-Golay method. 



Figure|7]illustrates the smoothing of an electrocardiogram signal by filtering the noise with Savitzky- 
Golay filter (third plot) and by denoising and enhancement with the fractal scalar conservation law ([B 
(second plot). In the second case, we can see that the relative maxima and minima are amplified. Of 
course, we can obtain even greater amplification of low/medium frequencies by suitably tuning the pa- 
rameters but, in these conditions, we will obtain sizeable variations between each peak. Thus, when we 
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Noiseless EGG vs. Savitzky-Golay filter 









500 1000 1500 2000 2500 3000 



Figure 7: Top: a noisy electrocardiogram u{x); Middle: the noiseless signal v{x) (blue) and denoising & 
enhancement of the ECG signal (red) using the FFT on the fractal equation ([T|) with An'^a = 0.005, b = 
0.0225 and A = 1.7. Bottom: filtering using Savitzky-Golay method. 



wish to amplify the low frequencies, we have to be careful that we do not alter too much the shape of the 
noiseless signal. 

To rate the performance of denoising and enhancement of our filter, we consider signals with different 
SNR, which corresponds to the power ratio between a signal and the background noise. 

In figure [H we took a small SNR. We can observe that shapes are amplified and that noise has been 
reduced significantly. But, we can also remark the undershoots just before and behind the shapes. This 
phenomenon has been highlighted in fl\ in the setting of dunes moiphodynamics. It is a consequence 
of the mass conservation property, see equation (fTOl) . We obtain similar results in figure |9] with a large 
SNR. Indeed, as we can see, the proposed method of filtering yields both an interesting denoising and an 
amplification of low/medium frequencies. The third plot conveys the filtering using the Savitzky-Golay 
approach. As agreed upon, we obtain a preservation of relative maxima. Moreover, comparing the output 
obtained thanks to our nonlocal method with the one of Savitzky-Golay filter, we notice the better ability 
to increase local extrema (contrast enhancement) while keeping a good denoising. This statement is 
confirmed by Figure [T0(b)[ where the signal has a medium SNR. We obtain a "good" smoothing and an 
interesting amplification of the low and medium frequencies. 



13 



Noisy signal vs. filtered signal 
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Noiseless signal vs. filtered signal 
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Figure 8: First plot (top) represents a noisy signal u{x) (blue). Second plot (bottom) illustrates the 
noiseless signal v{x) (blue) and denoising & enhancement of the signal (red) using the FFT on the 
fractal equation O with Air'^a = 0.1, b = 0.3 and A = 1.5. 



5 Performance metrics 



In this section, we wish to measure the denoising ability of our model. To evaluate our approach and 
compare it with the Savitzky-Golay filter, we use two measures: the Mean Squai^e EiTor (MSE) and the 
Signal-to-Noise Ratio (SNR). These metrics are frequently used in signal processing. They ai^e defined 
as follows: 



MSE 



1 



N 



i=l 



SNR 



.dh - lOlogio ( J-"' ^ 



where uq is the noiseless original signal, u is the filtered signal and N is the length of the filter. 
It is easy to see that a small MSE corresponds to a high noise reduction and that a large SNR indicates 
a good denoising. To compare the performance of filters, we consider a signal of trigonometric type and 
an ECG signal. Noise is added to these signals with SNR varying between to 8. Results are plotted in 
figures [T0(a)| and [5(b)[ For each signal, we use a sample of 100 random noises. The performance of the 
two approaches is estimated using MSE and SNR criteria. We plot an average of the results on figures 
nn and [121 Figures |1 l(a)| and |12(a)| show SNR values for these two methods applied to trigonometric 
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Noisy signal vs.filtered signal 
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Noiseless signal vs. Savitzky-Golay filter 
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Figure 9: Choice of parameters: 47r^a = 0.01, b = 0.05 and A = 1.5. The red line depicts the filtered 
signal. 
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Noiseless signal vs. filtered signal 
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Noiseless signal vs. Filtered signal 
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(a) A-Tv^a = 0.01, b = 0.03 and A = 1.5 



(b) Aiv^a = 0.01, b = 0.05 and A = 1.5 



Figure 10: The red line depicts the filtered signal. In this case, we took uo{x) = cos(57rx) + cos(207rx). 



and electrocardiogram signals. Figures [T 1 (b) | and 1 12(b) | show MSE values. These values show the high 
performance of nonlocal approach in signal denoising for trigonometric signals. Remember that in this 
case, the algorithm used for the implementation of the equation is the FFT, which is most suitable for 
trigonometric signals. We may just notice in figure [12] that for high SNR, the Savitzky-Golay filter is 
better than the nonlocal filter approach. Still, when the SNR is low - and it is the critical case - the 
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proposed method is more efficient than the Savitzky-Golay. Thus, the implementation of our PDE based 
on the FFT may not be effective for any type of signal, at least when the SNR is high. Alternatively, we 
may use the finite difference scheme. Results of the filtering are illustrated in figure [13] As we can see, 
the proposed model is always more efficient than the Savitzky-Golay approach, but for low SNR the FFT 
scheme remains best. Finally, regardless the type of signal considered, the fractal conservation law ^ is 
a good tool of denoising, provided that it is implemented with the right method. 




(a) SNR values in Savitzky-Golay method in comparison (b) MSB values in Savitzky-Golay method in comparison 
to proposed method to proposed method 

Figure 11: Choice of parameters: 47r^a = 0.01,6 = 0.03 and A = 1.5. The initial signal is uo(x) = 
cos(57r2;) + cos(207rx). 




(a) SNR values in Savitzky-Golay method in comparison (b) MSB values in Savitzky-Golay method in comparison 
to proposed method to proposed method 

Figure 12: Choice of parameters: Air'^a = 0.005, b = 0.015 and A = 1.7 . The initial signal no(x) is an 
electrocardiogram (ECG) signal. 

Remark 6. Let us briefly explain how SNR metrics allow us to optimize the choice of parameter b. 
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(a) SNR values in Savitzky-Golay method in comparison (b) MSE values in Savitzky-Golay metliod in comparison 
to proposed method using finite difference scheme to proposed method 

Figure 13: Choice of parameters: a = 0.5, b = 0.3, A = 1/3. The initial signal uo{x) is an electrocar- 
diogram (ECG) signal. 




(a) (b) 

Figure 14: Choice of parameters: An'^a = 0.01 and A = 1.5 . The initial signal is uo{x) = cos(57rx) + 
cos(207rx). 

In figure \T4l we display the behaviour of SNR values for different values of b. As we can remark, the 
denoising will be most efficient for b ~ 0.03. Hence visualization of the surface (figure \r4(a)\ and curves 
(figure \r4(b)\ enables us to find easily the parameter bfor optimal denoising. 

6 Concluding remarks 

Our first aim was to introduce a fractal conservation law for denoising and contrast enhancement of sig- 
nals. This device permits to reduce considerably the noise and to increase contrasts simultaneously. The 
study showed that our filter eliminates the high frequencies and amplifies the low/medium frequencies. 
In this paper, we also discussed the choice of pai^ameters a, b and A. 

This equation has been implemented using both a finite difference scheme and the fast Fourier trans- 
form. Various well-known measuring metrics have been used to compare our method with the Savitzky- 
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Golay filter. Results showed the good performance of our model. Moreover, the analysis highlighted 
that, depending on the considered signal, it may be more suitable to use the finite difference scheme or 
the FFT algorithm when the SNR is high. Obviously, for a sinusoidal type signal, it is preferable to use 
the FFT, whereas for a signal like step functions it is better to implement the fractal equation with finite 
difference method. Nevertheless, no matter the algorithm used, the fractal consersation law ([B is an 
interesting and natural method for denoising and contrast enhancement. 

These satisfying properties for signal processing encourages us to implement it for image enhance- 
ment: this study is in progress. 
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